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The main objective of spintronics is to un¬ 
derstand the mechanisms by which it is possi¬ 
ble to achieve efficient electrical control of spin 
configurations and of spin currents [H [2]. In 
the last decade, the way to achieve this ob¬ 
jective has experienced a breakthrough, due 
to: i) the discovery and understanding of a 

mechanism to generate spin currents in conduc¬ 
tors with magnetic order and in paramagnetic 
conductors/semiconductors [3ll5] . ii) the experi¬ 
mental observation of theoretically proposed spin 
injector systems [Tj [6l ^ and Hi) the synthe¬ 
sis of 2D materials with long spin relaxation 
timefTl [ 5 . The generation of spin currents, 
spin injections and spin conservation are medi¬ 
ated by the Spin-Orbit Coupling (SOC) mainly 
via topological effects and/or Rashba effect [9l 411| . 
Thus, the search for systems experiencing these 
properties is a primary concern for the develop¬ 
ment in spintronics. In this work, we propose 
a non-centrosymmetric honeycomb-lattice quan¬ 
tum spin hall effect family formed by atoms of 
the groups IV, V and VII of the periodic ta¬ 
ble, as shown in Fi^. We made a structural 
analysis and a Z 2 characterization. This material 
presents a giant Rashba-type spin-splitting with 
an unusual spin texture and a hexagonal warping 
effect, which lead to scattering process different 
from the usual materials with non-trivial topo¬ 
logical phases. We proposed a four band effective 
model to explain both the origins of such unusual 
spin texture in the Rashba effect and the SOC 
band inversion in the topological insulator phase. 
This dual behavior only has been reported in the 
BiTel 3D topological insulator |121 ITS] . 

Topological insulator materials support helical metal¬ 
lic edge states (or surface states for 3D systems), forming 
topological Dirac fermions protected by the time-reversal 
symmetry on an insulating bulk[3 [7]. This topological 
protection is responsible for the absence of backscatter- 
ing process as well as the in-plane chiral spin texture in 
this metallic states. This particular edge band disper¬ 
sion and spin helicity have been experimentally charac¬ 
terized via Angle-Resolved Photoemission Spectroscopy 
(ARPES) and localized Scanning Tunneling Microscope 
(STM) transport measurements [6l [Mj. On the other 
hand, the Rashba SOC leads to parallel spin-polarized 


band dispersion curves with in-plane opposite helical spin 
texture[T5]. allowing the control of the spin direction 
through an electric field [3]. This interaction originates 
from the lack of spatial inversion symmetry along the 
perpendicular direction to a two-dimensional electron gas 
(2DEG), which is hosted on the surface of a system or 
on an interface. The dispersion curves and Fermi con¬ 
tours have been characterized by spectroscopic measure¬ 
ments for many surfaces and interfaces [19]. Large 

Rashba couplings are found in materials formed by heavy 
elements with strong intrinsic SOC such as Bi, Pb, W and 
others fl^H^ . 

The intense study of the nontrivial topological phase 
in condensed matter systems began with the theoretical 
proposal of the quantum spin Hall effect in graDhene[24j. 
Despite its small SOC gap, this honeycomb-lattice ma¬ 
terial has promoted several searches for new materials 
with the same electronic and structural properties |25L 
1^ . Several graphene-like materials proposed theoreti¬ 
cally have been synthesized. However, most of these sys¬ 
tems, different from the quantum wells, have direct band 
gap in the K point of the Brillouin Zone (BZ). There¬ 
fore, the sub-lattice symmetry breaking bring the system 
to a trivial band gap, this effect can easily arise when the 
deposition on a substrate is considered. The functional¬ 
ization of germanene |33j and tin monolayers |28j , which 
exhibit a non-trivial band gap located in the F point, 
motivated us to search honeycomb-lattice systems with 
strong SOC and band gap located at this point. 

Geometry and electronic properties: We propose 
a system with a honeycomb lattice formed by atoms of 
the groups IV, V and VH of the periodic table, as shown 
in Fi#. Interpreting the hexagonal lattice as two tri¬ 
angular sub-lattices A and B, the system has a V atom 
type on the sub-lattice A, whereas in the sub-lattice B 
it has a IV-VH dimer. These non-centrosymmetry sys¬ 
tems have a buckled format and fulfill the symmetry op¬ 
erations of the C^y point group symmetry. From the 
ab initio calculations of the phonon spectrum for this 
family of materials, we find that the system formed by 
Bi, Pb and I atoms is mechanically stable. The forma¬ 
tion energy is about —0.77 eV, which is calculated as: 
Ef = EpbBii - Atpb - MBi - Ml: where EpbBii is the en¬ 
ergy of the system, and MPb,Bi,i are the chemical potential 
from the respective bulk phase. 

In Fig[T] we present the band structure and the pro¬ 
jections of the wave function in the Pz and Px^y atomic 
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Figure 1. Side and top view of the atomic structure of PbBil (a), and Band structure with SOC (b) and without SOC (c). The 
projection in the + Py orbitals of all atoms are indicated by the blue circles, while the Se p^ orbitals are indicated by the 
red circles. Both color are in gradient scale depending on weight of that projection in each state. 


orbitals for each band and fc-point, for both cases tak¬ 
ing into account and disregarding the SOC. First, we will 
discuss the band structure and the atomic orbitals inter¬ 
action disregarding the SOC. We find that the system 
is a direct gap semiconductor. The valence band maxi¬ 
mum and the conduction band minimum are located in 
the gamma point. Since the system has no magnetic mo¬ 
ment, the Time Reversal (TR) symmetry is preserved 
and the spin components are bi-degenerate, as shown in 
Fig[^. The valence orbitals of IV, V and VI atoms are 
mainly p-type orbitals, these orbitals are responsible for 
the chemical bonds. The crystal has a layered structure, 
so the z direction is different from the x or y directions 
in the atomic plane. Thus, there is an energy splitting 
between and Px^y orbitals. We find that the band 
structure near the Fermi energy is described by p-type 
Bi atomic orbitals. For Bi atom the Px,y orbitals have 
lower energy than pz orbitals. Consequently before SOC 
is considered, the conduction band mainly consists of Pz 
Bi orbitals, whereas the valence band is dominated by 
the Px^y Bi orbitals. 

Next we include SOC effect in the above atomic pic¬ 
ture. The atomic SOC Hamiltonian is given by Hso = 
XS ■ L with A = 2 mXc^r ^ depending on the detailed 
potential U of atoms, which couples orbital angular mo¬ 
mentum to spin [34). From the symmetry properties and 
how this Hamiltonian modify the atomic spectrum of Bi, 
Pb and I atoms in the crystal, we find that the spec¬ 
trum energy at the F point is given by the {lAj,^^)} 
effective states. J is the effective total angular moment 
with projection jz in the z axes and A corresponding to 
the Bi and Pb-I contributions. Transforming the Px and 
Py orbitals to p± with definite orbital angular momen¬ 
tum, we have J = 3/2,1/2 and jz = ±3/2, ±1/2. The 
|A3/2,±3/2) effective states are associated with the p_|_ 
and p_ effective orbitals, whereas the |Aj,±l/2) effec¬ 
tive states are linear combinations of p orbitals, in this 


case mainly of Pz orbitals. Thus, as shown in Fig[J, at 
the F point, the character of the band structure is in¬ 
verted, the conduction band mainly consists of Px,y Bi 
orbitals while the valence band is dominated by the Pz 
Bi orbitals. The inclusion of the SOC change the gap 
direction and the band structure is highly modified. The 
valence band maximum is located in the F —^ AT direction 
very close to the F point and the conduction band mini¬ 
mum is located in the M —>■ AT, as shown in Fig[T]:. We 
classify the SOC effects in the band structure into two 
classes: i) at the time reversal invariant momenta such 
as F and M the degeneration remains. Far from these 
points, the spin degeneracy is lifted and the M-shape as¬ 
sociated with the Rashba effect appears in the valence 
bands. This effect is characterized by the energy split¬ 
ting, Eh, represented in the Fig[^; ii) by turning the 
strength of the SOC, Asoc, where Asoc = 0(1) means 
absence (full presence) of SOC, we verify that at the F 
point a topological phase transition occurs. The system 
pass through a metal state with Asoc = 0.65, showing 
a band inversion characterized by energy gap at the F 
point, Egp, which it is not always equal to the band gap 
of the system Eg. Later, we will discuss both Rashba 
effect and topological properties in the PbBil system. 

In the lowest energy configuration, the Bi-Pb and Pb-I 
{h) bounds have 3.04 A and 1.35 A of length, respec¬ 
tively. Whereas, the buckled, d, have 1.3 A of length. 
Using GGA functional, we find that the lattice constant 
is 4.77 A, whereas for LDA is 4.53 A, as shown in Fig[|i. 
We find that the substitution of Pb atoms by Sn atoms or 
changing the lattice constant allows the control of both 
the band gap direction (Fig[^) and the like-Rashba spin¬ 
splitting in the valence band (Fig[^). The decrease (in¬ 
crease) of the lattice constant can be understood as a 
compressive (tensile) strain. We noticed three different 
behaviors depending on the intensity of the strain ap¬ 
plied: semi-metal, indirect band gap semiconductor and 
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Figure 2. (a) Energy in function of lattice parameter depending on the exchange-correlation functional used, i.e Local Density 
Approximation (LDA) or Generalized Gradient Approximation (GGA). (b) Gap size in function of lattice parameter, where 
Eg and Egr are the gap size of the system and only in the F point respectively. The dashed line at 0.0 of the gap represent the 
value where the variable Eg become inverted, (c) Size of Rashba spin splitting in meV in function of lattice parameter. 


direct band gap semiconductor, as represented in Fig[^ 
by the regions I, II and III, respectively. The equilibrium 
position is represented by arrows, shown as in the Fig[^ 
As mentioned above, in the region II, in which the system 
is in the equilibrium lattice parameter, Egr is negative 
indicating a band inversion at the F point. The gap Egr 
is different from the indirect band gap Eg , which was also 
considered negative for comparison purposes. In the re¬ 
gion I, when we apply compressive strain, the minimum 
of the conduction band, in M —>■ iF direction, begins to 
descend and cross the Fermi energy for a = 4.3 A, form¬ 
ing a semi-metal. To represent this behavior, we change 
the sign of the band gap Eg. In the region III, we consider 
a tensile strain. Eg and Egr are equals and the system 
behaves as a semiconductor with a direct band gap at the 
F point. For any value of strain, the states near to the 
Fermi energy in the F point remain inverted, so if this 
band inversion is associated with a topological state, this 
state will be robust against strain, even for a tension of 
20%(5.6 A ), as shown in Fig I 

Unconventional Rashba effect: According to our 
calculations, the value of the spin-splitting in the valence 
band due to the Rashba effect is about 60 meV, which is 
huge compared with the values found in 2D systems and 
surprisingly is on the order of the highest found in 3D 
systems |35j. As mentioned above, this phenomenon is 
robust to strain, even for large compressive strain, large 
Rashba values {Er = 90 meV) may be obtained, as 


shown in Fig [I- Although the conduction band does 
not have the M-shape characteristic of the Rashba ef¬ 
fect, the (Sz) spin texture evidence the same properties 
of the systems that exhibit this effect. Specifically, in 
both valence and conduction bands, the usual (Sz) spin 
texture associated with the Rashba spin-splitting is ob¬ 
served, as represented in Fig[^. However, the and 
(Sy) spin components exhibit an unusual behavior. In 
Fig [3)3, we show the in-plane spin texture in six con¬ 
stant energy planes. These planes cut the two conduc¬ 
tion bands named CBl and CB2, and the two valence 
bands named VBl and VB2 (see Fig §1). Near to the 
F point, the VBl exhibits a clockwise helicity, as shown 
in the energy plane el. Far from this point (|fe| greater 
than 0.1 in the reciprocal space), this helicity switches to 
counterclockwise. Very close to the F point, the VB2 ex¬ 
hibits a counterclockwise helicity. This helicity switches 
to counterclockwise for fc-values greater than 0.01, as is 
observed in the e2 and e3 energy planes. For BVl and 
BV2 the spin texture behaviors can be divided into three 
regions of the fc-space, |fe| being the norm of the vec¬ 
tor k = {kx,ky): i) when |fe| is less than 0.01, the spin 
texture of both bands rotates in the same direction, ii) 
when |A;| is greater than 0.01 and less than 0.1, the bands 
have the same helicity, unlike the usual behavior in the 
Rashba effect and in) when k is greater than 0.1, the 
bands have different helicity. The in-plane spin rotation 
in the CBl and CB2 have the same behavior described 
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Figure 3. (a) Bandstructure of PbBil in the direction K’ - F - K, with six cut el — 6 in constant energy of the valence and 
conduction bands VBl, VB2, CBl and CB2 respectively, (b) Spin texture from cuts of the bandstructure el —6. Spin projection 
in {x, y) plane are represented by the size of the arrown. For {Sz) we used a red-white-blue color gradient scale. In both figures 
dark red and dark blue represent a {Sz) up and down spin configuration respectively. 


above. However, the CBl and CB2 have opposite helicity 
direction to that found in the VBl and VB2, as shown 
in e4, e5 and e6 energy planes in Fig[^. 

We observed a hexagonal modulation of the bulk states 
analogous to the existence of a nonlinear term in the 
SOC[3S], as shown in the e6 and e2 energy planes in 
Fig[^. The nonlinear terms in the SOC lead to higher 
order components of the Rashba effect [37], which in¬ 
duce this hexagonal warping effect. This modulation is 
also observed in the surface states of the 3D topolog¬ 
ical insulators and is associated with the symmetry 
breaking [34l I36j . From our knowledge, there are no other 
2D systems experiencing this hexagonal warping effect 
in the bulk states. Additionally, in the spin texture, we 
observed a periodic {Sz) spin precession, which is asso¬ 
ciated with the lack of space inversion symmetry and 
the threefold rotation symmetry breaking. Interestingly, 
although the band structure describes a 2D semiconduc¬ 
tor, the {Sz) spin polarization increases as the energy 
increases. Thus, near to the Fermi energy the {Sz) spin- 
polarization is null, the same as occurs in the surface 
states of Bi2Te3[35]. The consequence of these {Sz) spin- 
component fluctuation and hexagonal warping effect in 
the electronic and transport properties has been exten¬ 
sively studied |36l [38l [39] . 

In the Rashba effect, at any energy plane, the eigen¬ 
state tpnik) has opposite spin direction to the eigenstate 
tpn{—k) with the same band index and opposite mo¬ 
menta. Thus, near to the Fermi energy, the backscat- 
tering processes are forbidden. This behavior is also ob¬ 
served in the PbBil system, as shown in Fig[^ On the 


other hand, in the Rashba effect, the inelastic scatter¬ 
ing between different bands is possible, i.e, the eigen¬ 
state ijjn{k) has the same spin direction to the eigenstate 
fe -|- Afc). Because of the particular in-plane spin 
texture helicity in the bulk states of the PbBil system, 
this scattering process in the PbBil system is different. 
As is observed in Fig Wp, except for the small energy 
range characterized by the e2 energy plane, in order for 
inelastic scattering processes between different band to 
occur, the {Sx) and {Sy) spin flip is necessary. Provided 
the TR symmetry is preserved, this protection against 
backscattering will remain. These ”TR protected bulk 
states” has no origin in the topological properties of the 
system and is only a result of the combination of the ef¬ 
fects caused by the high order Rashba components. This 
will be discussed later through the analysis of each of 
these components. As far as we know, the PbBil system 
is the first system that has such spin texture properties 
in their bulk band structure. This behavior is of deep 
interest in spintronics. 

Topological insulator phase: We hnd that in addi¬ 
tion to the inherent TR protected bulk states, the PbBil 
system is also a 2D topological insulator with ID TR pro¬ 
tected edge states. This topological protection against 
any TR invariant perturbation, can be seen through the 
topological order described by the Z 2 invariant miiD]- 
[42]. Thus, in order to fully determine the topological 
character of the non-centrosymmetric PbBil system, we 
implemented the evolution of Wannier Center of Charges 
(WCCs) as an alternative method for the Z 2 invariant 
calculation using ab-inito simulations (see methods). In 
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Figure 4. Spin texture band structure of armchair nanoribbon of the system PbBil with the {Sx)(a.) left, {Sy){a,) right, (S'z)(b) 
spin components projection. The width of the calculated ribbon is 98.7 Ain (c) we have the evolution of Wannier Center of 
Charges, and if we draw an horizontal line of reference that cross the picture, we would have at least an odd number of crossing, 
showing that this material is a topological insulator. 


Fig [^, we show that there is no a horizontal reference 
line that cross the evolution of the WCCs at least an odd 
number of times, showing that this material is a topolog¬ 
ical insulator [ 45 . 

Using electronic structure calculations, we verify that 
the band structure of a nanoribbon of the PbBil system 
exhibits edge states in the middle of the bulk band gap, 
regardless of the nanoribbon orientation. In Fig |^, 
we present the spin texture in the band structure of 
an Armchair nano-ribbon, which has two Dirac cones 
in the F point formed by bands with opposite spin, as 
expected. When the nano-ribbon is constructed, the 
mirror symmetry is broken and therefore, the Dirac 
cones located at different edges are nondegenerate. As 
shown in the Fig[^ although the edges states satisfy the 
TR symmetry and have the chiral spin helicity, the spin 
has preferentially the kx direction. The Sy and Sz spin 
components are lower for both bulk states and metallic 
edge states. We hnd that in the nanoribbon band struc¬ 
ture, the protected bulk states are maintained, but now 
these are unidirectional. Thus, in the nanoribbon, the 
backscattering is forbidden for both edge conductivity 
channels and bulk conductivity channels. This suggests 
that, contrary to what happens in most systems with 
nontrivial topological phases, the bulk states would not 
pose a problem for spintronic devices. 

Model: In order to understand both topological prop¬ 
erties and unconventional spin texture in the Rashba ef¬ 
fect, we propose a simple model that captures the elec¬ 
tronic properties near to the Fermi energy. Because the 
valence and conduction bands are formed mainly by the 
effective states {|Bi,/, j^)} with J = 3/2,1/2, we use these 


states to form the basis of our model. By construction, 
this base contains intrinsically the full SOC. From a phe¬ 
nomenological analysis of the electronic structure, we ini¬ 
tially describe the effects, in principle, required to repro¬ 
duce the results obtained via DFT calculations. Subse¬ 
quently, using symmetry elements, we construct a tight- 
binding effective model, showing that these effects arise 
naturally when the fc —>■ F limit is considered. 

Although, the Rashba effect and the surface states in 
3D topological insulators arise as a result of different 
physical properties, the effective Hamiltonian that cap¬ 
tures the first order Rashba effect in 2D semiconductors 
systems is also used to describe the surface states in a 3D 
topological insulatorj^: 'H(fc) = /lo(fc) + hRx{k). The 
first term /lo(fc) = is the quasi-free electron approxi¬ 
mation and the second term (fe) = or-^ [cr x k) ■ z rep¬ 
resents the interaction between the spin Pauli matrices 
vector cr and the momentum k. Here or-^ is the Rashba 
coefhcient, and z the unit vector perpendicular to the 
plane in which lies the semiconductor. When the term 

which depends on the effective mass m* and the elec¬ 
tron mass rrie sufficiently small relative 

to parameter or, states with linear dispersion and he¬ 
lical spin texture appears, describing the surface states 
in a 3D TI. Otherwise, the Rashba effect band disper¬ 
sion appears. Since |l/2, ±1/2) effective states are those 
that form the valence band, in which the M-shape of the 
Rashba effect occurs, only these states have the hR^{k) 
term. Obviously both the |l/2,±l/2) as the |3/2, ±3/2) 
effective states have the ho{k) term, which introduces 
the band parabolic effect. The hexagonal warping on all 
bands is introduced considering the three order Rashba 
term hR^{k) = aR^[{kx + iky)^ + {kx-iky)^]az^\- Join- 
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Figure 5. In (a) we have the band structure from the DFT approach (blue lines) and the effective tigh binding model (circles). 
In (b) we calculated the spin component extracted from the model. The values for each variable of the model is given in the 
table see methods). 


ing these terms, the effective Hamiltonian can be written 
as: 


\ nl,{k) n,/,{k) J ’ 


( 1 ) 


where T-Li/ 2 {k) = hy^{k) + h]l^{k) + h](^{k) and 

T~l- 3 / 2 {k) = —hy^{k) + h^j^^{k) are the effective Hamilto¬ 
nian that described the 11/2, ±1/2) and 13/2, ±3/2) effec¬ 
tive states, respectively. The phenomenological analysis 
of the behavior of the band structure does not provide in¬ 
formation on the interaction Hamiltonian between these 
states, 'Hint{k). Therefore, we propose a fuller analy¬ 
sis through a tight-binding effective model. The tight- 
binding Hamiltonian matrix elements in the momentum 
space is written as. 
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H{k) = e{k)+ ^ (2) 

(0=i» 


where tn- = {n, J,is the nearest neighbor 
hopping term, with n indicating the lattice site and 
corresponds to the i-th of the six nearest neighbor vector. 
Since the |(Pb-I)j, jij) effective states are not considered, 
we omit the terms associated with the near neighbors. 
Based on the symmetry properties of the system, the 
hopping terms can be uniquely determined. We study 
the behavior of the tight-binding Hamiltonian matrix el¬ 
ements in the fc —>■ T limit, obtaining the most simplified 
model. From this model we have the relations between 
the ttijj and parameters and the hopping terms 
and the Hint interaction term (see methods), explicitly: 
^ = 3a^Re(t^^), = |lm(t^^), = ^a^Im(t^^) 

and Hint = icrzjaRe{k^){kj; +iky). 


The Hint interaction is critical to the spin polariza¬ 
tion transfer and the induction of the helical spin tex¬ 
ture for the band where the major contribution comes 
from the |i?i3/2,±3/2) effective states. Thus, the effects 
mentioned above together with the band inversion allow 
the reproduction of both band structure and spin texture 
obtained via DFT calculations, as shown in Fig[^ The 
procedure followed for this Hamiltonian is described in 
the methods and the optimized values were obtained via 
a least squares standard approach in order to match DFT 
calculation (See Table [I]). 


TTli 7Tl2 Oiint 


1/2 

3/2 

Values 0.3714 0.5294 3.5853i 

3.0919 

12 

20 


Table I. Values of the parameters for the effective tigh binding 
model. 


CONCLUSIONS 

From ab initio calculations we made a crystalline struc¬ 
ture analysis, a Z 2 topological evaluation and a mapping 
of topological edge states. The Rashba spin splitting re¬ 
ported here is up to 60 meV without strain, and rise 
up to 90 meV with strain, providing to this system a gi¬ 
ant Rashba effect. We also noticed an exotic spin texture 
with spin flipping, an unconventional effect in the Rashba 
splitting where, at constant energy planes in the same 
band index, we have two circular spin helicity turning in 
opposite ways, an effect of interest in spintronic devices. 
We also concluded that PbBil is a topological insulator 
via the approach of the switching of Wanniers Centers 
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of Charges. This material is strong upon strain, main¬ 
taining the topological properties even for a tension of 
20 %, although compression may leave the material in a 
semi metal electronic phase. We mapped the topological 
edge states with opposite spin as expected from TRS, not 
perfectly degenerate at the inversion point, i.e T, due to 
the lack of inversion symmetry. Our finding may enable 
further research in this new topological material, which 
present such unusual spin texture. 


METHODS 

DFT calculations; The first-principles calculations 
were performed within the Density Functional Theory 
(DFT) framework, as implemented in the SIESTA 
code|44). considering the on-site approximation for the 
spin-orbit coupling|45j. The Local Density Approxi¬ 
mation (LDA)[1^ is used for the exchange-correlation 
functional. We used an energy cutoff of 410 Ry to dehne 
the grid in the real space, and the supercell approxima¬ 
tion with a k-points sampling for the reciprocal space 
integration equivalent to 20 x 20 x 1 in the unitary 
cell. The 2D PbBil sheets lie in the xy plane, and a 
vacuum of 20A was used in the z-direction to avoid the 
undesirable interaction between the periodic images of 
PbBil sheets. The system geometry was fully relaxed 
until the residual forces on the atoms were smaller than 
0.01 eV/A. 

Z 2 invariant calculation: The easiest way of dis¬ 
covering new topological materials is the evaluation of 
an invariant number always associated with a particu¬ 
lar symmetry. For systems protected by Time Reversal 
(TR) symmetry, the invariant characterizing the system 
is called a Z 2 invariant. When this number is not zero, 
the system has TR protected metallic edge states. 

The formulation of the Z 2 number follows many differ¬ 
ent formulations: via the parity eigenvalues for systems 
with inversion symmetry 133, the formulation via Resta 
polarization |42j . by an integration of the Berry phase 
over the effective Brillouin zone [ID], and by the switch¬ 
ing of Wannier Center of Charges (WCCs) [13] .When the 
materials are non-centrosymmetric, the ab initio calcula¬ 
tions can give a generic gauge [43] on the states. There¬ 
fore, we chose to calculate the Z 2 invariant via the switch¬ 
ing of WCCs. In that manner, we know that the Z 2 
invariant can be written as: 

Z2 = Pe(i = T/2)-Pe(i = 0), (3) 

where t represents an adiabatic time. Now we want to 
rewrite equation @ in terms of the WCCs. For that, 
we need to define Wannier functions (WF). The WFFs 
belonging the the same unit cell R have the following 


format: 

|Pn) = (4) 

A WCC Xn is defined as the mean value of (0n|A|0n) 
where X is the position operator and |0n) is the state 
corresponding to a WF in the cell with i? = 0. So we 
have: 

i r 

= — y dk{unk\dk\unk) ■ (5) 

Assuming that we got ^ §bz ^ ~ 

/, //, where sumation in a represents occupied states and 
A the Berry connection, we finally have: 

^2 = E - ^"{TI2)\ - E [;si(0) - 5:i'(0)] . 

Ct Oc 

( 6 ) 

With the smooth gauge condition between t G [0,T/2] 
it is possible to track the evolution of WCCs during 
a half cycle, which is also smooth. It is important to 
note that this argumentation is only valid when some 
properties are valid for the construction of the Wannier 
Functions. Now if we have an even or odd number of 
crossing between half cycle of the WCCs, we are dealing 
with a trivial or a topological material, i.e Z 2 equal 0 or 1. 
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